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Abstract. We study one-parameter expanding evolution families of simply connected domains 
in the complex plane described by infinite systems of evolution parameters. These evolution 
parameters in some cases admit Hamiltonian formulation and lead to integrable systems. One 
example of such parameters is complex moments for the Laplacian growth that form a Whitham- 
Toda integrable hierarchy. Another example we deal with is related to expanding coefficient 
bodies for conformal maps given by Lowner subordination chains. The coefficients bodies are 
proved to form a Liouville partially integrable Hamiltonian system for each fixed index and the 
first integrals are obtained. We also discuss the contact structure of this system. 



1. Introduction 

Complex dynamical systems, iterations and construction of Lie semigroups with respect to 
the composition operation lead to the study of expanding systems of plane domains. These 
systems can be thought of as one-parameter families of domains (typically simply connected 
and univalent) that form subordination chains in the Riemann sphere C, which are defined 
for <t <T where T may be oo. This means that 0,{t) C ^{s), i}{t) ^ ^(s), whenever t < s. 

Kufarev ^S] studied a one-parameter family of domains ^l{t), and regular functions g{z,t), 
g : Q{t) X [0, T) — >■ U, where U = {( : |C| < !}■ He proved differentiability of g{z, t) with respect 
to t for z from the Caratheodory kernel $7(0) of 17 (t). 

Sometimes, the evolution of ^}{t) can be characterized by infinite sets of evolution parameters 
given as a result of integration of certain dynamical systems. Hamiltonian interpretation of 
such systems yields their integrability. The interest to integrable systems is caused first of 
all by Hamiltonian systems of mathematical physics where the Hamiltonian is thought of as 
an energy functional and its critical points correspond to minimal energy of particle systems. 
However, other areas of mathematics also use Hamiltonian formalism, e.g., optimal control, 
classical mechanics, etc. 

One of the examples of integrable systems of evolution parameters corresponding to one- 
parameter families of expanding domains is the problem of Laplacian growth (Hele-Shaw prob- 
lem). The Laplacian growth means the Dirichlet problem for a harmonic potential where the 
boundary of the phase domain 0,{t) is unknown a-priori (free boundary), and in fact, is defined 
by the normality of its motion. The classical (strong) solvability of this problem suggests that 
phase domains 0(t) are bounded by analytic curves. It turns out, that given an initial domain 
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il(0) the set of Richardson's moments (see the definition in |21]) solves completely the inverse 
potential problem, thus describing the evolution of phase domains as long as the classical so- 
lution exists. As it has been shown in jj^j, ^S|) HHi the Laplacian growth problem can 
be embedded into a larger hierarchy of domain variations (Whitham-Toda hierarchy) for which 
all the complex moments are treated as independent variables (generalized time variables), and 
form an integrable system. 

Another example of evolution parameters is the set of coefficients of the Riemann maps of U 
onto the family of the phase domains 0(t) that has attracted attention of the complex analysts 
for a more than eighty year period of the last century. The Bieberbach conjecture ||8j has proved 
to be the most intriguing problem that forced investigations in geometric function theory during 
this period. It finally has been solved in 1984 by de Branges [S], who proved that the n-th 
Taylor coefficient bn of a conformal homeomorphism /(C) , normalized by /(O) = 0, /'(O) = 1, 
does not exceed its index: < n. This class of functions we denote by S. The inequality is 
sharp and the equality is attained only for rotations of the Koebe function k{z) = z(l — z)~'^. 
De Branges' proof has put an end to the numerous attempts to prove or disprove the Bieberbach 
conjecture. However, a much more general and difficult problem is to describe the range of n 
first coefficients as a set of values = {(62, . . . , bn), / G S"}. In the trivial case V2 is the disk of 
radius 2. Only the first non-trivial coefficient body V3 = (b2,b3) has been described completely 
by Schaeffer and Spenser in 1951 in their famous monograph |^. A qualitative description of 
Vn, n > 3, has been partially given in [3]. Apart from these two monographs there are only few 
works where a progress in such a problem has been made (see, e.g., [23] )• 

The key stages of our presentation are as follows. Firstly, a boundary point of Vn is given by 
a unique function whose first coefficients b2, ■ ■ ■ ,bn define the rest of coefficients. Secondly, this 
function unlikely the Laplacian growth, presents a full mapping (the complement of the image of 
U is of measure zero) of the unit disk onto the complex plane minus an analytic graph with the 
root at infinity. Thirdly, erasing this graph we obtain a family of expanding domains completely 
described by the initial domain, and in its turn by the point of dVn- Finally, we obtain a (2n — 3) 
dimensional family of evolution dynamics completely described as a flow generated by a certain 
vector field, which admits a Hamiltonian formulation. It is proved to be partially integrable 
in the Liouville sense and the first integrals are obtained. Moreover, it turns out that the first 
integrals that generate additional directions form a horizontal space so this Hamiltonian system 
possesses a contact structure. It is important that the evolution exists all the time < t < 00. 

Our paper is devoted to a detail realization of this plan and the description of this system 
of evolution parameters, its Hamiltonian formulation, and finally, its integration. The central 
idea of the Hamiltonian formulation is forced by the Lowner evolution (see, e.g., |21 117| I2()j ) 
that is based on evolution equations in partial and ordinary derivatives mutually linked. The 
controllable dynamical system for evolution parameters turns out to be Hamiltonian under 
the assumption that the necessary optimality conditions hold. So, we are especially interested 
in subordination chains of domains generated by Riemann maps corresponding to the critical 
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controls, and therefore, to the trajectories, that in particular, give all boundary points to the 
coefficient body for univalent maps. 

From one side our main results describe evolution coefficient systems for certain type of 
Lowner's dynamics of domains and maps that admit the Hamiltonian formulation and are inte- 
grable. The coefficient bodies will form a certain type of hierarchy (dependent on the coefficient 
indices) and the first integrals will be derived. From the other side, we give all boundary points 
of the coefficient bodies Vn- 

2. The Lowner and Lowner-Kufarev evolution 

There are many papers and monographs devoted to the Lowner and Lowner-Kufarev equa- 
tions. However we revisit this theory from the point of view of correspondence between Lowner 
equations in partial and ordinary derivatives and univalence of solutions to the partial derivative 
Lowner equation. 

Let us consider a subordination chain of simply connected hyperbolic domains ^}{t) in the 
complex plane C, which is defined for < t < T, where T may be infinity. We suppose that 
the origin is an interior point of the kernel of {Q{t)}J^Q. Let us normalize the growth of the 
evolution of this subordination chain by the conformal radius of ^}{t) with respect to the origin 
to be e*. By the Riemann Mapping Theorem we construct a subordination chain of mappings 
f{C,t), C ^U, where each function f{C,t) = e*C + t2(i)C^ + - • • is a holomorphic univalent map of 
U onto 0,(t) for every fixed t. Pommerenke |2n|l21j introduced such chains in order to generalize 
Lowner's equation. His result says that given a subordination chain of domains ^(t) defined for 
t G [0, T), there exists an analytic regular function 

p(C, t) = i+p^ (t)c + P2(t)C' + • • • , C(^u, 

such that Re p{C, t) > and 

(1) ^^ = ^^^^^^'*)' 

for Q GU and for almost all t € [0, T). The class of holomorphic functions p defined as above 
we denote by C. The equation is called the Lowner-Kufarev equation due to two seminal 
papers by Lowner with 

iu{t) I /- 

(2) P^^^t) = -^^, 

where u{t) is a continuous function regarding to t G [0,7"), and by Kufarev [TB| in general case, 
where this equation appeared for the first time. 

The equation represents a growing evolution of simply connected domains. Let us consider 
the reverse process. Given an initial domain Q(0) = VIq (and therefore, the initial mapping 
/(C,0) = foiO), and a function p(C)i) from the class C, we solve the equation (P) and ask 
whether the solution f{(^, t) represents a subordination chain of simply connected domains. The 
initial condition /(C,0) = /o(C) is not given on the characteristics of the partial derivative 
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equation hence the solution exists and is unique. Assuming s as a parameter along the 
characteristics we have 

dt dQ df 

as as as 

with the initial conditions t(0) = 0, ({0) = z, /(C) 0) = /o(C)) where z is in U. This leads us 
to the Cauchy problem for the Lowner-Kufarev equation in ordinary derivatives for a function 

C = w{z,t) 

dw , , 

(3) — = -wp[w,t), 

with the initial condition w{z, 0) = z. 

We see that the equation Q is exactly the characteristic equation for Unfortunately, 
this approach requires the extension of foiW^iC^t)) iiito U because the solution to is the 
function /(C, t) given as fo{w~^ {(, t)), where C = w{z, s) is a solution of the initial value problem 
for the characteristic equation Q that maps U into U . Therefore, the solution of the initial 
value problem for the equation may be non-univalent. 

Let S stand for the usual class of all univalent holomorphic functions /(z) = z + a2Z^ + . . . 
in the unit disk. Solutions to the equation Q are regular univalent functions wi^z^ 
02(^)2^ + . . . in the unit disk that map U into itself. Conversely, every function from the class 
S can be represented by the limit 

(4) f{z)= lime'w{z,t), 

t— >oo 

where there exists a function p{z,t) from the class C for almost all t > 0, such that w{z,t) is a 
solution to the eqiation Q (see [23 pages 159-163]). Each function p{z,t) generates a unique 
function from the class S. The reciprocal statement is not true. In general, a function f € S 
can be determined by different functions p (z C. 

Prom |211 page 163] it follows that we can guarantee the univalence of the solutions to the 
Lowner-Kufarev equation in partial derivatives assuming the initial condition /o(C) given by 
the limit (jlj) with the function p{-,t) chosen to be the same in the equations and ©. 

In his original paper JJj Lowner dealt with functions that map U onto domains each of which 
was obtained by slitting C along a unique Jordan curve. This led him to a subclass S' of the 
class S of one-slit maps which is dense in S. Each function / E S" is obtained as the limit (HJ 
where w{z,t) is a solution to the Cauchy problem for the Lowner equation 

, . dw g»«{t) _|_ y; 

^ ^ 'dt ~ ~^e*"(*) -w' 

with the initial condition w{z, 0) = z where u{t) is a continuous control function. The corre- 
sponding Lowner equation in partial derivatives is just equation ^ with the function p{(, t) 
given by The solution to this equation is univalent all the time if and only if the 
initial condition is given as the limit for the same p given by d^J). It maps the unit 
disk onto a growing family of one-slit domains. The class of such domains lies dense in the space 
of all simply connected domains in the sense of Caratheodory's kernel convergence. 



UNIVALENT FUNCTIONS AND INTEGRABLE SYSTEMS 



5 



We proceed with generalized Lowner chains of domains with several slits. Let il.{t), € ^{t), 
t € [0, T), be a subordination chain of simply connected domains each of which is obtained 
by slitting C along a finite number of Jordan curves with m + 1 distinct endpoints. The m + 
1-th endpoint is at infinity. The one-parameter family f{C,t) stands for the corresponding 
subordination chain of normalized conformal maps f{C,t) = e*C + &2(i)C^ + - ■ ■ , C ^ U, t € [0,r). 
These maps satisfy the equation (0) with the function p{(, t) given by Q where u{t) is a 
piecewise continuous control function regarding to t G [0,T) (see, e.g., Goluzin 12 ). However, 
this representation is not unique and another approach has been suggested originally by Kufarev 
and described in It was based on the function p{C,t) in given by the following formula 



where Uk{t) are continuous control functions regarding to t € [0, T) and Afc(t) are measurable 
non-negative control functions, X]fcLi'^fc(*) = 1- 

Constructing the characteristic equation Q we come to the Lowner-Kufarev equation in 
ordinary derivatives and the limit (jlj gives a subclass S" of the class S of multi-slit maps which 
is also dense in S. The function w{z, t) in (jlJ) is a solution to the Cauchy problem for the 
Lowner-Kufarev equation 



with the initial condition w{z,Q) = z where tifc(t) and \k{t) are defined as in @. 

Similarly to the one-slit evolution, the equation with the function p given by © possesses 
univalent solutions all the time if and only if the initial condition is given as the limit Q) for 
the same p in ((TJ given by 

Although the representation (^IHl) or Q is not unique neither, it allows to elaborate a unique 
one. Given a subordination chain n[t) of multi-slit domains there exists a unique set of positive 
numbers {Xk}t=v YJk=i^k = 1, and a unique set of continuous control functions {ukit)}^=i^ 
such that Q.{t) = f{U,t), where / is a solution to the equation with the function p given by 
(|ni), and Afc(t) = A^. This statement follows from the results of An analogous uniqueness 
statement for the representation is true too. 



As it was stated in Introduction, one of the examples of integrable systems of evolution 
parameters corresponding to one-parameter families of domains is the problem of Laplacian 
growth. In its classical formulation it deals with the smooth family {Q{t)}, t E [0, T) of phase 
domains, i.e., when d^}{t) are smooth (C°°) interfaces for each t, and the normal velocity Vn at the 
boundary continuously depends on t at any point of dO,{t). Each 0,{t) is supposed to be simply 
connected in C for any t € [0, T) fixed. Starting from the initial moment t = the boundary 
becomes even analytic smooth for t S (0, T). Richardson's complex moments completely describe 
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this evolution, say given an initial phase domain Q{0) there exists a unique infinite set of the 
moments for Q{t) at each t G (0, T) and vice versa. One of the reasons for a slit dynamics to 
appear in this process is extending the solution to exist beyond the blowup time T. One way 
of the regularization of the problem is to introduce a surface tension parameter instead of zero 
boundary condition in the classical Hele-Shaw formulation (see, e.g., |14|'). This can produce a 
crack morphology in which a part of the boundary develops a thin finger penetrating the phase 
domain. Another proposition (see J3|) is to continue smooth dynamics by slit dynamics beyond 
the blowup time when a possible cusp is developed at the boundary. Of course, in this case 
one relaxes smoothness of the boundary permitting the corresponding Riemann map to have 
singularities at the unit circle. 

Another reason to consider slit dynamics is the determinate evolution, duality, and univa- 
lence of the solutions to the Lowner and Lowner-Kufarev equations (in ordinary and partial 
derivatives) stated in the preceding section. 

In view of the involved Riemann maps the most natural evolution parameters are the Taylor 
coefficients of these maps. Of course, the infinite set of coefficients completely describes the 
evolution given by the maps f : U ^ ^l{t). However in general, the differential equations for 
these coefficients obtained from the Lowner or Lowner-Kufarev equations do not admit any 
reasonable Hamiltonian formulation that leads to integration of this (infinitely dimensional) 
system. 

Instead we propose to construct hierarchies defined by the coefficient bodies as an object to 
study. Let / S 5. We construct the coefficient bodies by setting 

oo 

Vn = {{b2,...,bn) : f{0=C + J2^kC' &S}, n>2, 

k=2 

in the (2n — 2)-dimensional Euclidean space. Every point of the boundary of Vn is given by a 
unique function / G that maps the unit disk onto the plane minus piecewise analytic Jordan 
arcs forming a tree with a root at infinity having at most n — 1 tips. Each Vn defines a class {i^o} 
of initial domains for multi-slit subordination evolutions {Q(t)}, and the corresponding systems 
of differential equations for coefficients form an integrable system Z)„ for each n, moreover, its 
dual system L>* obtained by the characteristic equation admits Hamiltonian formulation. The 
latter systems form hierarchies {-Dn}^2 {■^n}^=2- This idea up to some extent will be 
proved in the next sections. 

4. Coefficient bodies 

By the coefficient problem for univalent functions we mean the problem of precise finding 
the regions Vn defined above. These compact sets have been investigated by a great number of 
authors, but the most remarkable source is a famous monograph written by Schaeffer and 
Spencer in 1950. Among other contributions to the coefficient problem we distinct a monograph 
by Babenko that contains a good collection of qualitative results on the coefficient bodies Vn- 
The results concerning the structure and properties of Vn include 
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(iii) 



IV 



(i) Vn is homeomorphic to a (2n— 2)-dimensional ball and its boundary dVn is homeomorphic 
to a (2n — 3)-dimensional sphere; 

(ii) every point x € dVn corresponds to exactly one function f G S which will be called a 
boundary function for Vn ; 

with the exception for a set of smaller dimension, at every point x € dVn there exists a 
normal vector satisfying the Lipschitz condition; 

there exists a connected open set Xi on dVn, such that the boundary dVn is an analytic 
hypersurface at every point of Xi. The points of dVn corresponding to the functions 
that give the extremum to a linear functional belong to the closure of Xi. 

It is worth to note that all boundary functions have a similar structure. They map the unit 
disk U onto the complex plane C minus piecewise analytic Jordan arcs forming a tree with a root 
at infinity and having at most n — 1 tips. This assertion underlines the importance of multi-slit 
maps in the coefficient problem for univalent functions. 

Lowner's approach is based on the following idea: a function mapping U onto the complement 
of a single slit admits univalent dynamics for t E [0, oo) represented by the Lowner equation (^ISJ 
forming a subordination chain f{C,t) = X]JS=i ^n(^)C"- Let us deduce a system of differential 
equations for the coefficients of f{C,t) substituting the expansion of f{C,t) into (0). This gives 



kbi 



fc-i 



-i{k-j)u 



bk{0) = ak, t>0, 



e*. 



A; = 1,2, ... , where /(C,0) = /o(C) = E^i^nC"- In particular, bi{t) 

Going to the Lowner equation in characteristics © we write w{z,t) = Yl'^=i'^n{t)C^ and 
substitute this function in ((S)). To formulate the result we introduce the following notations 



/ aiit) \ 



ait) 



Ait) 



/o 

aiit) 
a2{t) 



Then the differential equation for a(i) is of the form 

n-l 

(8) d = -a - 2 ^ e"' 






. 





. 


aiit) 


. 


an-2it) ■ 


. aiit) 









a(0) = a° 



In particular, ai(t) = e * and limt^oo e*afc(t) 



/ 1 \ 
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5. HAMILTONIAN FORMULATION FOR THE COEFFICIENT SYSTEM. InTEGRABILITY 

5.1. Hamiltonian dynamics and integrability. Let us recall briefly the Hamiltonian and 
symplectic definitions and concepts that will be used in the sequel. There exists a vast amount 
of modern literature dedicated to different approaches to and definitions of intagrable systems 
(see, e.g., 0, (H, 0, EHl). 

The classical definition of a completely integrable system in the sense of Liouville applies 
to a Hamiltonian system. If we can find independent conserved integrals which are pairwise 
involutory (vanishing Poisson bracket), this system is completely integrable (see e.g., [1], 
[7j). That is each first integral allows us to reduce the order of the system not just by one, but 
by two. We formulate this definition in a slightly adopted form as follows. 

A dynamical system in C^" is called Hamiltonian if it is of the form 

(9) X = VsH{x), 

where denotes the symplectic gradient given by 

„ , d d d d 

V, 



r\— 5***5 r\— 5 r\ 5 * * * ) r\ 

dx„+i C)X2n OXi dXr. 

The function if in Q is called the Hamiltonian function of the system. It is convenient to 
redefine the coordinates (xn+i, . . . , x^n) = (tpi, ■ ■ ■ , tpn), and rewrite the system Q as 

OH — dH 
(10) Xfc = ^1;,^ = - — , k=l,2...,n. 

dijjk OXk 

The system has n degrees of freedom. The two-form uj = Yl^=i dx A dip admits the standard 
Poisson bracket {•, ■} 



associated with uj. The symplectic pair (€2*^,0;) defines the Poisson manifold (C^", {•, •}). These 
notations may be generalized for a symplectic manifold and a Hamiltonian dynamical system 
on it. 

The system H10|) may be rewritten as 

(11) Xk = {xk,H}, 'iljk = {ipk^H}, fc = l,2...,n, 
and the first integrals ^ of the system are characterized by 

(12) {^,H} = Q. 

In particular, {H, H} = 0, and the Hamiltonian function H is an integral of the system ©. If the 
system (jllj) has n functionally independent integrals $i, . . . , which are pairwise involutory 
^j} = 0, k, j = 1, . . . , n, then it is called completely integrable in the sense of Liouville. The 
function H is included in the set of the first integrals. The classical theorem of Liouville and 
Arnold 3- gives a complete description of the motion generated by the completely integrable 
system (|11() . It states that such a system admits action-angle coordinates around a connected 
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regular compact invariant manifold. One can work with real Poisson manifolds instead, making 
use of the real Haniiltonian function 2Re H keeping all other formulas changeless. 

If the Hamiltonian system admits only 1 < k < n independent involutory integrals, then 
it is called partially integrable. The case k = 1 is known as the Poincare-Lyapunov theorem 
which states that a periodic orbit of an autonomous Hamiltonian system can be included in a 
one-parameter family of such orbits under a nondegeneracy assumption. A bridge between these 
two extremal cases k = 1 and k = n has been proposed by Nekhoroshev JHI and proved later 
in ^Uj, The result states the existence of ^-parameter families of tori under suitable 

nondegeneracy conditions. 

5.2. Coefficient system. We will show that the system Q becomes an integrable system when 
treated as a description of the boundary hypersurface dVn- Solutions a{t) to (jS)) for different 
control function u{t) (piecewise continuous, in general) being multiplied by the factor e* represent 
all points of dVn as t ^ oo. The trajectories e^a{t), Q <t < oo, fill Vn so that every point of Vn 
belongs to a certain trajectory e^a{t). The endpoints of these trajectories can be interior or else 
boundary points of V^. In this way, we set Vn as the closure of the reachable set for the control 
system ©. 

According to property (ii) of Vn given in the previous section, every point x € dVn is attained 
by exactly one trajectory e^a{t) which is determined by a choice of the piecewise continuous 
control function u{t). The function f & S corresponding to x is a multi-slit map of U. If the 
boundary tree of / has only one tip, then there is a unique continuous control function u{t) 
in t G [0, oo) that corresponds to /. The case of multi-slit maps we will consider in the next 
section. 

To reach a boundary point x G dVn corresponding to a one-slit map, the trajectory e^a{t) has 
to obey extremal properties, i.e., to be an optimal trajectory. The continuous control function 
u{t) must be optimal, and hence, it satisfies a necessary condition of optimality. The Maximum 
Principle is a powerful tool to be used that provides a joint interpretation of two classical 
necessary variational conditions: the Euler equations and the Weierstrass inequalities (see, e.g.. 



EH)- 



To realize the maximum principle we consider an adjoint vector 



V Mt) J 



with complex valued coordinates ipi, . . . ,ipn, and the real Hamiltonian function 
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where ip means the vector with complex conjugate coordinates. To come to the Hamiltonian 
formulation for the coefficient system we require that ^ satisfies the adjoint system of differential 
equations 

Taking into account Q we rewrite as 



n-l 



(14) ^=i^E + 2Y,e—^Hs + l){A'rji^, 

where E is the unit matrix. 

The maximum principle states that any optimal control function u*{t) possesses a maximizing 
property for the Hamiltonian function along the corresponding trajectory, i.e., 

(15) maicH{a*{t),^*{t),u) = H{a*{t),ib*{t),u*), t>0, 

u 

where a* and ip* are solutions to the system (|H1 114p with u = u*{t). 
The maximum principle (jlSf) yields that 

dH{a*{t),ip*{t),u) 



(16) 



= 0. 

u=u* (t) 



du 

Evidently, ©, (0, and ^ imply that 

dH{a*{t),rit),u*it)) _ 
^ ' dt 
for an optimal differentiable control function u*{t). 

The Hamiltonian formalism for system ((5l ll4() will lead to integrability. First, we show how ip 
can be expressed in terms of the phase variable a. In the following theorem causing no confusion 
we will use *(...) to denote the matrix transposition. 

Theorem 1. Let a{t) and ■0(t), ipiT) = ^{vi,... ,f„) obey the system |3 \T^, T > 0. Then 
i^kit) = Cn-k+i, k = 1, . . . ,n, where ci, . . . ,Cn are the Taylor coefficients of the expansion 

{vnz + ■■■ + viz'')w'{z,T) ^ 

— — 

Proof. Let w{z,t) = Xlfc^i (^k{t)z'^ be a solution to the Lowner differential equation Q. Differ- 
entiating © with respect to z we immediately have 

d f z \ z fe'^' + w le'^'w 

(1^) 31 ,..//-. - -JmK -17, — - + 



dt \w'{z,t) J 'w'{z,t) \e™ — w (e™ — w)^ 
Considering the expansion 
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we obtain 
(19) 



dqjt) 
dt 



n-l 



s=l 



where q{t) = . . . ,qn{t))'^. We observe that system (fT^ differs from system (fTl|) only by 

the transposition sign. 

In order to satisfy the condition q{T) = {vi, . . . ,Vn)'^ we denote by 

(v„zH + viz'')w'{z,T) ^ ^ 

si^^ t) = = g ckit)^ > 

and observe that g{z,t) obeys the same equation (fTK|) where z/w'{z,t) is substituted by g{z,t). 
Hence, 

c{t) = 

obeys the system (fT^ substituting q{t) by c{t). It is easily seen that 

/ Vr,. \ 



c{T) 



The difference in the transposition sign implies that ^ipk = Cn-fc+i, k = 1, . . . , n. This completes 
the proof. □ 

Putting T = in this theorem we come to the following corollary. 

Corollary 1. Let a{t) and ^^{t), V'(O) = (vi, . . . , Un)^ obey the system t > 0. Then 

i'kit) = Cn-k+1, k = 1, . . . ,n, where ci, . . . ,Cn are the Taylor coefficients of the expansion 



{VnZ H h viz"') 

w'{z,t) 



Y.^k{t)z\ 



k=l 



Remark 1. Since Im ai{t) = 0, the Hamiltonian function H{a, -0, u) does not depend on Im ipi. 
Therefore, without loss of generality, we can put "01(0) = vi to be real. This assumption leave a, 
■02, • • • , "iAn changeless. 

Now we are able to apply the maximum principle. An optimal continuous control function 
satisfies the maximizing property (|15|) and obeys equation (|16|) . The Hamiltonian H{a,'ip,u) is 
a trigonometric polynomial with respect to u of degree n — 1 if 7^ 0. Let "0(0) = (^i) • • • , VnY' ■ 
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Note that = e*. We assume that v„ 7^ 0, and diminish nton — 1, otherwise. At t = 0, we 

have 

n-1 

H{a^, {vi,. . . , VnY" , u) = — ■uie"* — 2 y~](Re Vg+i cos{su) + Im ^5+1 sin(su)). 

s=l 

For the optimal u*, 

n-1 

Hu{a^, {vi, . . . , Vn)'^, u*) = 2 ^^(s Re Vg+i sin(s'u*) — s Im Vg+i cos(s'u*)) = 0. 

s=l 

Suppose that V2, ■ ■ ■ ,Vn are such that 

n-1 

Huu{a^, (ui, . . . , Vn)'^, u*) = 2 ^(s^ Re u^+i cos(sit*) + Im Vg+i sin(su*)) 7^ 0, 

s=l 

for the optimal u*. As we see, the term with vi does not influence the process of maximizing, 
and we will assume it zero later in course. The inequality holds for all {v2, ■ ■ ■ ,Vn) G M^""^ 
except for a set A of dimension at most 2n — 4. Indeed, due to the maximizing property of the 
optimal control function, this condition breaks down if 

Huu{aP, {vi,..., Vnf, u*) = Huuu{a°, {vi,---, Vnf , u*) = 0. 
Thus, the set A is determined by two linearly independent equations 

n-1 

Re Vs+i cos{su*) + Im Vs+i sin(s'u*)) = 0, 

s=l 
n-1 

y~^(s^ Re Vs+i sin(s'u*) — Im Vg+i cos(s'u*)) = 0, 

s=l 

with the fixed optimal u*, varying V2, ■ ■ ■ ,Vn- Therefore, A is a linear space of dimension 2n — 4 
in M2n-2_ 

Definition 1. We say that a vector {v2, ■ ■ ■ ,Vn) satisfies the regularity condition on [0,T] if 

Huu{a*{t),r{t),u*)^0, te[0,T], 

for the optimal u* and = (ui, . . . j^n)"^ along the optimal trajectory {e^ a* (t) , tp* (t)) corre- 

sponding to u* . 

Evidently, if {v2, ■ ■ ■ , Vn) G M^"~^ \ A, then (^2, . . . , Vn) satisfies the regularity condition on 
[0, T] for T > and small enough. Let us denote by ^(r) the set of all (f2, • • • satisfying 
the regularity condition on [0,r]. 

Let {v2,...,Vn) e Y{T). Then 

(20) Huu{a,^,u)j^O, te[0,T], 

for {a,ijj,u) from a neighborhood of {a* , i/:* , u*) and this neighborhood depends on t. The 
equation 

Hu{a,il;,u) = 0, 
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together with the regularity condition on [0, T] determines an analytic implicit function u = 
u{a,ilj) in a neighborhood of {a*,'tp*), and u{a,ip) lies in a neighborhood of u* = u{a*,ip*). 
We substitute u = u{a, ■0) in the phase system Q and the adjoint system ((T^ and solve 



the Cauchy problem with the initial data a(0) = a , ip{0) = {vi, 



Since u = u* along 



(a*,'i/'*), the Hamiltonian H[a,ip,u) satisfies the maximum principle. Note that the neigh- 
borhood of {a*{t),ij)*{t)) can be determined by a neighborhood of the corresponding vector 
{v2-, ■ ■ ■ iVn) = (V'2(0), • • • , V'n(O)). According to Corollary Q the function il){t) can be rep- 
resented as a function of the phase variable a and of the initial conditions say 
V'(i) = Lp{a{t),vi, . . .,Vn). 

Theorem 2. Let (t>2, . . . , Vn) satisfy the regularity condition on [0, T]. Then the system (0 \14\ ) 
with u = u{a,ip) = u{a, (p{a, 0,V2, ■ ■ ■ ,Vn)) is partially Liouville integrable. Moreover, the first 
integrals form a contact structure. 

Proof. We substitute ip{t) = ip{a{t), vi, . . . , Vn) in the Hamiltonian and obtain 

H{a{t),'ip{t),u) = H{a{t),ip{a{t),vi, . . .,Vn),u) = n{a{t),vi, . . .,Vn,u). 

To each point (vi, . . . ,Vn) from a neighborhood of a vector satisfying the regularity condition, 
there corresponds an optimal continuous control function. The maximum principle implies that 



TCu{a{t),vi, ...,Vn,u) 



0, 



u=u{a,ifi{a,vi,...,Vn)) 



that, together with ((T7|) . gives 

(21) H{a{t),vi, . . . ,Vn,u{a,ip{a,vi, . . . ,Vn))) = const. 

In order to prove the partial integrability we will find the first complex integrals ^"1, . . . $n- 
We already have one real integral (|2H) . 

By the corollary from Theorem ^ we deduce that 



^k+iz^ + w'{z,t) Ckz''. 



(22) '^Vn-k+iz'' = w'{z,t)Yi>„ 

k=l k=l k=n+l 

We denote by *('^i, . . . , ^*n) the vector of the first integrals of the Hamiltonian system (jSl [Tl|l 
given by 

/ ^1 \ / 

ai 2a2 ... (n - l)a„_i nan 

^•2 ai . . . (n - 2)a„_2 (n - l)a„_i 

= ... (n-3)a„_3 (n - 2)a„_2 



■02 



J \ 0... ai J yiPn J 

with the control u given in the statement of the theorem. Indeed, the equality (|22|) implies that 
^fc = ''^fc are constants for all t and k = 1, . . . ,n. The first integral H21|) allows us to conclude 
that 7i} = 0. One checks that {^s, ^k} = -s^s+k-i for ah 1 < s < A; and s + k - 1 < n. 
Otherwise, $fc} = for s + k — 1 > n. This implies that 
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• [n/2] first integrals {*&[n/2]+i; • • • > *^n} are pairwise involutory; 

• the integrals {$i, . . . , *l*[n/2]} ^.re not pairwise involutory but their Poisson brackets give 
all the rest of integrals. This structure is said to be contact (which is not uncoupled 
because, for example, {$1, . . . , $„} = —^n 7^ 0). 

This completes the proof. □ 

Remark 2. We observe that 

n-l 

H{aQ, (ui, . . . , Vn)^ , u*) = —vi — 2 ^^(Re Vs+i cos(su*) + Im u^+i sin(su*)) 

s=l 

att = 0, which is the constant in the right-hand side of \21\) . It contains an additive term (— fi). 
The optimal control function u* does not depend on vi . Hence, we can put ui = without loss of 
generality. This assumption does not influence a,'ip2, ■ ■ ■ md defines uniquely the constant 
in the right-hand side of i21\) . Thus, the first integral \21\) can he rewritten as 

(23) H{a{t),V2, . . . ,Vn,u{a,<p{a,V2, ■ ■ ■ ,Vn))) = const. 

The integral H23|) gives a local parametric description of the boundary of the coefficient body 
for bounded univalent functions from S (such that f £ S, \f{z)\ < M, I < M < e'^). The 
points of this part of the boundary correspond to the functions that map U onto the disk 
\w\ < M slit along an analytic curve with exactly one tip at its interior point. The boundary 
hypersurface is homeomorphic to an open set on the 2n— 3-dimensional sphere in M^"~^ and must 
be parameterized by 2n — 3 free parameters. We should select 2n — 3 independent parameters 
among 2n — 2 variables V2, ■ ■ ■ ,Vn- Observe that the Hamiltonian H{a,ip,u) and the solution 
^ to the adjoint system ()13() are linear with respect to ip. The implicit function u = u{a, ip) is 
invariant upon multiplying ■0 by a positive number. Therefore, both H and ijj may be determined 
up to a positive multiplier without loss of generality. Putting, e.g., \vn\ = 1 we reduce the number 
of free parameters in H23|) to 2n — 3 that gives a local parametric representation of the boundary. 

The integral (|21|) is continued from [0, T] to a bigger interval [0, T+e) as long as the regularity 
conditions are preserved. Let T > be the minimal positive number such that the regularity 
conditions break up at T iov V2 ■ . . ,Vn, i.e., 

'Huuia(T),V2, . . . ,Vn,u{a{T),ip{a{T),V2, . . . ,f„))) = 0, 

'Huuuia(T),V2, . . . ,Vn,u{a{T),ip{a{T),V2, . . .,Vn))) = 0. 

Denote by 

Gl(t,f2, ... ,Vn) ■■= 'Huu{a{t),V2, . . . ,Vn,u{a{t),if{a{t),V2,. . . ,Wn))), 
G2(t,f2, ...,Vn):= 'Huuu{a{t),V2, . . . ,Vn,u{a{t),ip{a{t),V2, . . . ,Vn))). 

In these formulas we admit that a{t) is determined by t>2, • • • , fn according to the phase system 
0. The system Gi{t,V2, . . . ,Vn) = G2{t,V2, ■ ■ ■ ,Vn) = of linearly independent equations 
defines a (2n — 4)-dimensional manifold in the (2n — 3)-dimensional space of free parameters 

V2, ...,V„, \V„\ = 1. 
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Thus, all integrals (^5]) are continued on t £ [0, oo) for all V2, - ■ ■ ,Vn except for a (2n — 4)- 
dimensional set. Using continuity of integrals (|2,S|) with respect to V2, ■ ■ ■ ,Vn we conclude that 
the relation ^ holds for all iv2,...,Vn) € M^^-s ^nd for ah t > 0. 

6. MULTI-SLIT EVOLUTION 

Now we generalize the results of the preceding section to the multi-slit evolution. We do this 
for the maps with the boundary tree having exactly two finite tips because other cases can be 
easily obtained by analogy. The Lowner-Kufarev equation has the form 



dw 



+ W 



—— = —W X— 

dt \ e*"H*) - w 



+ (1-A)- 



+ w 



iu2(t) 



W 



where ui and U2 are continuous control functions and A G (0, 1). The differential equation for 
phase variables ((H} becomes 

n-l 



-ISUl 



+ {1- X)e-''''^)A'a, 



n-l 



(24) d = -a-2^(A 
The Hamiltonian is written as 

H{a, ip, ui,U2, A) = Re 
and the adjoint system H14() as 

(25) ^ = Ie + 2 ^(Ae-^^"i + (1 - A)e-^^"2)(s + 1)(^'^)' ) ^. 



2 ^(Ae-*^"i + (1 - X)e-''''^)A'a 
s=l J 



n-l 



s=l 



The maximizing condition H15() now splits into two bits 
(26) max H{a*{t),^*{t),ui,U2,X) = H{a*{t),i;*(t),ul,u*2,X), t>0, 

ui,U2,X 

where a* and ip* are solutions to the system (1241 125p with ui = u1{t), U2 = U2it), and moreover, 



(27) 



Re 



0. 



Indeed, if condition (|27|) were not true the value of the optimal A would be or 1, and the 
corresponding boundary function would be a one-slit map that contradicts our supposition. The 
maximum principle (|26|) implies that 



(28) 



dHia*it),rit),ui,U2,X) 



0, j = l,2. 



Uj=U*{t) 
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Evidently, (El, (ESI), and (UHl) imply that 

for optimal differentiable control functions ul{t) and U2{t). 

Obviously, Theorem ^ and the corollary thereafter still true in this case and tp can be ex- 
pressed in terms of the phase variable a. Continuous optimal control functions u^, satisfy 
the maximizing property (|2())) and obey equations (|28|) . The Hamiltonian H{a,tlj,ui,U2, \) is 
a trigonometric polynomial with respect to ui and U2 of degree n — 1 if ipn 7^ 0. Let again 
-0(0) = {vi, . . . , u„)^, and we assume that u„ 7^ 0, and diminish n to n — 1, otherwise. We have 
V'n(i) = e*. An additional parameter is A € (0, 1), which is constant. At t = 0, we have 

Hia^, (fi, . . . ,Vn)'^ ,ui,U2, A) = Aif(a°, {vi, . . . ,Vnf ,ui) + (1 - \)H{a^, {vi, . . . ,f„)^,U2). 

The vector of parameters is (A,fi, . . . ,Wn)j and we put f 1 = by the same reason as in the 
preceding section. Let us consider (A, f2, . . . , w„) G M^"^-*^ (or better (0,1) x M^"~^). For the 
optimal u\ and U2 

Fui(a°, (fi, . . . ,Vnf ,u\,u*2, A) = Hu2{a^, (^i, • ■ • ,Vnf ,ul,u*2,\) = 0. 

Suppose that A, f2, . . . , are such that 

{vi, ■ ■ ■ ,Vnf ,u\,u*2,\), or Hu^u^ia^, {vi, . . . ,Vn)'^ ,ul,U2, X), 

do not vanish. This holds for all {X,V2, ■ ■ ■ ,Vn) G C (0, 1) x M^""^ except for a set A of 
dimension at most 2n — 4. The linear space B is obtained by restriction ()27|) which, being 
rewritten at t = for the optimal control functions, is equivalent to 

i/(a°, {vi,..., vnf, ul) - H{a\ (t^i, . . . , t;„)^, ul) = 0. 

So, the dimension of 5 is 2n — 2. The set A is determined by two linearly independent equations 

HuiuiicL^, {Vi,... ,Vn)^ ,ul,U2,\) = Huimuiia'^ , (vi,.-. ,Vn)'^ ,ul,U2,\) = 0, 

or 

HuiUiia'^, (Wl, . . . ,Vn)'^,ul,U2, A) = Hu2U2U2{a^, {vi, . . . ,Vnf ,u\,u*2, A) = 0, 

with the fixed optimal U2, varying [\,V2-, ■ ■ ■ iVn) in B. Therefore, ^ is a linear space of 
dimension 2n — 4 in S C M^""^. 

Definition 2. We say that a vector (A, U2, . . . , Vn) satisfies the regularity condition on [0, T] for 
multi-slit evolution if 

-ff«i«i(a*,'0*,'"i,'"2i'^) / 0> or Hu2U2{a*,'>P*,u\,u*2,\)^Q, tG[0,r], 

for the optimal u\, U2 and V'*(0) = (wi, . . . , Un)"^ along the optimal trajectory (e*a*(t), V'*(i)) 
corresponding to u\, U2- 
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If (A, ^2, . . . ,Vrt) £ B \ A, then (A, ^2, . . . ,Vn) satisfies the regularity condition on [0,T] for 
T > and small enough. Let us denote by Y(T) the set of all (A,?;2,... ,Vn) satisfying the 
regularity condition on [0, T]. 

Let {v2, • • • , Vn) € Y{T). Then 

Huiui{a,tp,ui,U2,X) 0, or ^„2„2(a, ■0, t^i, ^2, A) / 0, t £ [0,T], 

for (a, tp, ui,U2, A) from a neighborhood of {a* , ul,U2, A) and this neighborhood depends on 
t. The equations 

Hu-^{a,ip,ui,U2,X) = Hu2{a,ip,ui,U2,\) = 0, 
together with the regularity condition on [0, T] determine analytic implicit functions ui = 
ui{a,tp,X) U2 = U2(a, V'jA) in a neighborhood of (a*,'0*,A), and each of ui = ui{a,tp,X) 
U2 = U2(a, "0, A) lies in a neighborhood of u\ = ui{a*,ip* , X), U2 = U2{a, tp, A) respectively. 

We substitute ui = iii(a, V', A), U2 = U2{a,ijj,X) in the phase system and the adjoint 
system and solve the Cauchy problem with the initial data o(0) = a*^, V'(O) = (f 1, • • • , Vn)"^ , 
{X,V2, ■ ■ ■ , Vn) € B. According to Corollary Q the function tp(t) can be represented as a function 
of the phase variable a and the initial conditions {vi, . . . , say ip{t) = if{a{t),vi, . . . , Vn, A). 
As in the preceding section, both H and ■0 may be determined up to a positive multiplier without 
loss of generality. Putting, e.g., = 1 we reduce the number of free parameters to 2n — 3 that 
gives a local parametric representation of the boundary. Repeating all further steps of preceding 
section we prove the following theorem. 

Theorem 3. Let (A, f2, . . . ,Vn) satisfy the regularity condition on [0, T]. Then the system \24\) 
with 

ui = ui{a,ip, X) = ui{a,ip{a{t),vi, . . . ,w„. A), A), 

U2 = ni(a, -0, A) = U2{a,ip{a{t),vi,. . . , w„. A), A), 

is partially Liouville integrable. Moreover, this statement is continued for all {X,V2, ■ ■ ■ ,Vn) G B 
and for all t > 0. 

Theorems 2,3 also solve the problem of integrability for the corresponding coefficient systems 
in Section 4 substituting the optimal control functions in the Cauchy problem. One can find 
exact solutions in explicit form for n = 3 in, e.g., |5] and j26j . 
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